Bifurcation between 2-component and 3 -component 
ground states of spin-1 Bose-Einstein condensates in 

uniform magnetic fields* 

Liren Lin^'^ and I-Liang Chern^'^'^'^ 

en ■ ^Department of Mathematics, National Taiwan University, Taipei 106, Taiwan 

, Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, Taiwan 

^SJ ' ^ Center of Mathematical Modeling and Scientific Computing, Hsinchu 300, Taiwan 

D 

^. 
T— I ! 

""I ' Abstract 

We study antiferromagnetic spin-1 Bose-Einstein condensates under external unifomi 
^ i magnetic fields. The bifurcation between two-component and three-component regimes 

I of ground states with respect to the magnetization M and the quadratic Zeeman effect 

q is justified. The proof is based on the technique of "mass redistribution" introduced 
in the authors previous work, which gives interesting inequalities and equalities satisfied 
by ground states. Some open problems arising naturally from our investigation are also 
• discussed in the end. 

a^ 

P 1 Introduction 

' Ever since the first realization of Bose-Einstein condensation (BEC) in 1995 lUl IH |9l], 
it has drawn great attentions of physicists as well as mathematicians. In early experi- 
ments, the atoms were confined in magnetic traps, in which the spin degrees of freedom 
^ ■ are frozen. By the mean-field approximation, such a system is then described by a 
■ scalar wave function, which satisfies the Gross-Pitaevskii (GP) equation JH [121 . In 
contrast, in an optically trapped atomic BEC, all its hyperfine spin states can be active 
simultaneously, and a spin-F BEC has to be described by a (2F + 1) -component vector 
function \E' = {ipF,ipF_i, ■ ■ ■ .ip^pY . Such spinor BEC was first realized in a gas of 
spin-1 23Na atoms in 1998 by the MIT group [1251 |26l |22l 12 El, and soon after that its 



theory was developed independently by several researchers [[231 [131 [IS. Since then, it 
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has also become a subject of intensive studies, both theoretically and numerically. On 
the other hand, although there are already many mathematical investigations of BEC sys- 
tems with multiple (even general A^) components [12T|. the specific structures of spinor 
BEC are not paid full attentions to by mathematicians. In this paper we consider the 
mean field ground states of spin-1 BEC. Before introducing the problems and goals to 
be studied, we shall describe the mathematical model first. 

As mentioned above, in the mean field approximation a spin-1 BEC is described by 
a three-component complex- valued ^ = (^i, V'o, ^-i)^- Since we will only be interested 
in ground states, we consider ^ as a function of the space variable x g M"^ and is inde- 
pendent of time. Under a uniform magnetic field, the system is described by the energy 



functional mm 



where V{x) is a real- valued function, (3n,(3s,p,q are real constants, and F = {Fx, Fy,Fz) 
is the triple of spin-1 Pauli matrices, which are given by 




F, = 




( 1 










Thus ^!*F^! denotes the vector {^!*Fx^, '^*Fy'^, '^*Fz'^). The notation |^| denotes the 

Euclidean length , and similarly for |V^j| and |\E'*F\I/|. We remark that 

here the coefficients of the terms of E[^] are normalized and are different from those 
used in the literature for simplicity. Physically, V{x) represents a state-independent trap 
potential, the terms with coefficients and f3s describe the collisions of the atoms, 
and p, q give the linear and quadratic Zeeman effects, which can be tuned by changing 
the magnitude of the applied magnetic field. The number of atoms A/'[\E'] and the total 
magnetization M of the system are given by 



1^1 



l^-il^). 



And a ground state is a minimizer of E under fixed M and M. Thus it's a varia- 
tional problem with two constraints. By normalization, we shall assume Ul^] = 1, 
and M[^] = M for some constant M. Note that |A^[^]| < Af[^] for every state \&, so 
we must have \M\ < 1. Due to the symmetry of the roles of ipi and ip-i, we will only 
consider < M < 1. We will consider the following setting for our model, which is 
quite general for the interested phenomenon in this work. 
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(Al) V e Lf^^, and V{x) tends to infinity unifonnly as x tends to infinity. Precisely 

lim (essinf = oo. 

R^oo \x\>R 

In particular V is bounded from below. 
(A2) f3n> (repulsive) and /9s > (antiferromagnetic). 
(A3) g > 0. 

Assumption (Al) will guarantee that V{x) traps the repulsive system in a localized re- 
gion, which is essential for the existence of ground states. Besides these three assump- 
tions, also note that due to the conservations of M and M, ground states are not changed 
by shifting the values of V{x) and p by any constants, and hence we shall also assume 
for simplicity 

(A4) y > andp = 0. 

This work is mainly motivated by the following observation: For fixed < M < 1, as 
q increases from zero, the ground state ^ undergoes a bifurcation from ipo = to ipo 
at a critical point qc{M), hence from a two-component (2C) profile to a three-component 
(3C) one. This phenomenon has been known for many years from numerical simulations 
Il28l [T9ll (A clear diagram showing the 2C regime and the 3C regime with respect to M 
and q is provided by Fig. 5 in [fTQl .). and was recently observed experimentally [fT4|. 
However, there seems to be no rigorous mathematical justification so far. In theoretical 
explanations by physicists, some simplified assumptions such as uniformity (assuming 
\1/ is a constant vector) or single-mode approximation (assuming the components of ^ 
are proportional to each other) are made, which have no sufficient reason for being 
satisfactory demonstrations. In [|20l, the authors found a principle which says that after 
a redistribution of the masses between different components, the kinetic energy will 
decrease. Using this fact, we successfully proved, among other things, the vanishing of 
ipQ of the ground state at g = (i.e. no external magnetic field). In the present work, we 
will show that the bifurcation phenomenon can also be deduced from the same principle, 
while not as obviously as before. Some basic properties about ground states have to be 
established first. It's interesting that many of the facts addressed in this paper can also 
be derived by using the idea of mass redistribution. 

The outline of the paper is as follows. Section 2 contains the basics of the model. We 
first introduce a well-known reduction by which one can simply study the amplitudes 
of the components of \E'. Then, after defining some notations in Section 2.1, the most 
fundamental properties such as existence, regularity and maximum principle are given 
in Section 2.2. In Section 3 we recap the idea of mass redistribution and, with the aid of 
it, prove some more useful facts about ground states: the continuity and monotonicity of 
ground-state energy with respect to M and g, the fact > \ip-i\, and the exponential 
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decaying of ground states as |a;| oo. In Section 4 we introduce an idea of perturbation 
by redistribution, which gives rise to some inequalities and equalities satisfied by ground 
states. In Section 5.1 the bifurcation phenomenon is justified by using the inequalities 
obtained in Section 4. Some approximations and characterizations of the bifurcation 
point induced by our proof are given in Section 5.2. In Section 6 we discuss three open 
problems naturally arising from this work. 



2 Preliminaries 



2.1 Reductions and notations 

We use B to denote the function space {H^(R^) n LliR^) n L'^(M3))3^ where Ll is the 



y-weighted space: A measurable function / is in L 



if 



/|/|V<oo. 



(When the domain of an integration is not specified, it's understood to be E"^.) B is a 
Banach space endowed with the following norm: 



1^1 



Given \& = {ipi, i/jq, ip^i) e B. Let ipj = \ipj\e^^^ be the polar form of the j-th component. 
Then we have (remember we have assumed p = 0) 



+ iV-^ip + 2|V>i||V>_i| cos(^i - 2^0 + O-i] 



+ (3s{\i^i\' - |V^-i|2)2 + + \dx. 



Since > |V|^j|| by the convexity inequality for gradients (see e.g. lUTl . 7.8), we 

find > S[{\^il l^ol, |^-i|)], where 



Slu] :-- 



for nonnegative triples u 
and satisfy 



(ui, Mo, u^i) e B. The equality holds if the 9j's are constants 



cos(^i - 2^0 + ^-i) = -1- 



(2.1) 
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Also, note that the conservations of M and M are constraints on and have nothing 
to do with the phases. These observations lead us to consider the following variational 
problem: 

Minimizing £ over the family of nonnegative triples 
u G B satisfying 7V[u] = 1 and A^[u] = M. 

It's not hard to prove that if ^ is a ground state, u = (l^/'il, |V^o|, iV'-il) is a solution of 
(*). Conversely, if u = (mi,mo,m_i) is a solution of (*), then ^ = V^o, V'-i) defined 



by ipj = uje^^^ for any choice of constants 9j satisfying (|2.1|) is a ground state. As 
a consequence, to understand the bifurcation phenomenon introduced in Section 1, it 
suffices to study (*). We will do so in the rest of this paper and no longer consider E. To 
facilitate later discussion we introduce some notations in the following. 

First let's make the rule that when a boldface letter, possibly with a superscript, is 
used to denote an element in B, its components are denoted by the same letter in normal 
font with indices 1, 0, —1. For example u denotes (mi, uq, as before, and similarly v 
denotes (fi,fo,f_i), denotes {w\,Wq,w^i), etc. We say a sequence in B converges 
weakly if it converges weakly in {H'^{^^)Y, in {L\^{^^)f, and in (L''(M3))3. Due to the 
reflexivity of H^, Ly and L^, every bounded sequence in B has a weakly convergent 
subsequence. Let 

B+ = {u G B I Uj > for each j } . 
We have the following fact. 
Lemma 2.1. B+ is a weakly closed subset ofM. 

Proof. Let {u'^} be a sequence in B+ which weakly converges to some g B. We 
need to prove > for each j. Let H+ = {u g (//^(E^))^ I > for each j }. Note 
that {u'^} is also a sequence in weakly converging to u°° in (i7^(]R3))'^. Since is 
a convex and closed subset of (if^(E'^))"^, if + is a weakly closed subset of (if^(E'^))3 by 
Mazur's theorem (see e.g. flSj, Theorem 3.7). Thus u°° g that is > for each 
J. □ 

The admissible class on which we are going to minimize £ is 
A = {u G B+ I M[u] = 1 and A^[u] = M} . 
The ground- state energy will be denoted by Eg, that is 

Eq = inf ^[vl. 

^ VGA 

And the set of minimizers is 

G= {u G A I ^[u] = Eg}. 
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For convenience, we will also frequently refer to elements in G as ground states. When 
M and q are regarded as variable parameters, we shall write Am, Eg{M, q), and GM,q to 
specify their values explicitly. Throughout this paper other parameters of our model are 
considered to be fixed and satisfy (Al) - (A4). Finally, we will use H{\i) to denote the 
integrand of £^[u], i.e. i^[u] = j H{\\). For convenience we write 

H{u) = HMn{^) + HpotiVi) + /fn(u) + Hs{u) + ff^ee(u), 

where Hkini'^) = |Vup, Hpot{u) = V|up, Hniu) and Hs{u) represent the terms with 
coefficients f3n and f3s respectively, and if^ee(u) represents the term with coefficient q. 
Accordingly, we also use Skin, ^pot, etc. to denote the corresponding parts of £. 

2.2 Basic properties 

In many aspects our three-component system can be regarded as a generalization of the 
one-component BEC model studied in IfTSl (though the main interests in this paper are 
different from theirs). The most fundamental properties about the one-component model 
(Theorem 2.1 of [flSl ) hold naturally for our three-component system (with however a 
remarkable exception: uniqueness. See Remark 12.21 and Section 6.1). We summarize 
them in the following theorem. 

Theorem 2.2. G 7^ 0, i.e. there does exist a ground state. Let u g G, then u is at least 
once continuously differentiable, and u satisfies the following Euler-Lagrange system in 
the sense of distribution: 

(/X + \)ui = Cui + 2/3s [uq{ui - u-i) + ui{u\ - + qui (2.2a) 
/iMo = Cuq + 2/3sMo(mi — u-if' (2.2b) 
(/X - A)m_i = + 2/3s \uq{u^i - ui) + - u\)\ + gw-i, (2.2c) 

where C = —A + V + 2/3^1 up, and X, /i are Lagrange multipliers arising from the con- 
straints ^f[u] = 1 and M[u] = M respectively. Furthermore, for each uj, either Uj = 
or Uj > on all ofM.^. 

We give some remarks. First, due to the repulsive assumption (A2), the existence 
result can be proved by the standard direct method in the calculus of variations, in which 
one tries to show that a minimizing sequence has a subsequence that weakly converges 
to an element in G. The only difference from a typical situation is that here the system 
is on the whole space but not a bounded domain. As a result, we do not have compact 
embedding to guarantee that the weak limit is still in the same admissible 

class. Instead, we should use the assumption (Al), which implies that most part of 
u is really contained in bounded domains, on which compact embedding applies. A 
precise argument can be given almost the same as in Lemma A.2 of IfTSil (see also [El|6l). 
Nevertheless, besides the conclusion of existence, some observations from the proof will 
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also be needed later. We give them in Lemma |231 below. The most important point is that 
we actually have strong convergence but not only weak convergence for the subsequence 
of the minimizing sequence. This holds for our model since all the terms (namely i^fcm, 
Hpot, etc.) of H are nonnegative. For convenience we give the proof in the appendix. 

Lemma 2.3. Let {u"} be a sequence in 1B+. Suppose A/'[u'*] — 1, A^[u^] — )■ M, and 
£[\i^] is uniformly bounded in n, then {u"^} has a subsequence {u'^^^^}^^ converging 
weakly to some u°° g A, which satisfies £^[u°°] < liminffc^oo ^^[u'^^'^^]- If we assume 
further that £[vi^] — )■ Eg, then G G, and u"(^) — )■ u°° in the norm ofM. 

Next, for the remaining assertions of Theorem 12.21 The Euler-Lagrange system (2.2) 
is called a time-independent Gross-Pitaevskii system (GP system for short). We remark 
that by definition (2.2) holds in the sense of distribution, while by approximation we 
see it's also valid when tested by elements in B. In fact, £,Af,M are continuously 
differentiable functions on B, and (2.2) is exactly 



where ' denotes Frechet derivative. We omit the routine verification of this fact. Once 
(2.2) is obtained, that u g G is continuously differentiable follows standard regularity 
theorem (e.g. [UTl, 10.2). And the strict positivity of a nonvanishing component can be 
obtained by the strong maximum principle (see e.g. [fTOl . Theorem 8.19). We will use 
these two facts tacitly to avoid repeatedly referring to Theorem 12.21 

Corollary 2.4. LetueG.IfO<M<l, then Uj ^ Q for j = 1,-1. 

Proof. Since j {u\-v?l^) = M > 0,uij^ 0, and hence ui > 0. To prove m_i 7^ 0, assume 
otherwise, then (2.2c) gives u^ui = 0, and so uq = 0. Thus among the three components 
only Ml 7^ 0, which implies M = 1, contradicting to our assumption. □ 

Remark 2.1. The assumption < M < 1 is necessary. If M = 1, it's obvious that only 
Ml 7^ 0. For M = 0, see Proposition 13 .71 and the remark following it. 

As is mentioned in the introduction, we will investigate whether uq vanishes or not, 
as a property depending on the values of M and q. Let's here use to denote the class 
of all u G B+ such that uq = 0, and let A*™" = A n B^f °. Note that for u g A*^° the 
constraints are 



^Af'[u] + -M'[u] = -£'[u] 



(2.3) 




and 




Then we define 
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Obviously, if u g G is such that no = 0, then u e 

The assertions in Theorem 12.21 for G (existence, regularity, and positivity of a non- 
vanishing component) also hold for G^^'^. The Euler-Lagrange system for u e just 
consists of (2.2a) and (2.2c) with uq = 0. A particular feature of the two-component 
system is the following convexity property. 

Lemma 2.5. Given u,v g Mf". Let w g be defined by = {uj + v])/2for 
j = 1, -1, then £[w] < {£[u] + ^[v])/2. 

Proof. Let D = {£[u\ + S[v])/2 - £[w], then D = Dkin + Dn + Ds, where 
which is nonnegative by the convexity inequality for gradients. Also, 

Dn = ^n[w] = — / (|u| -|v| ) >0, 



Ds = ^,[w] = — / (mi -fi >0, 



and 



as are easily checked. Thus > 0, which is what we want to show. □ 

From this convexity property we obtain the following uniqueness result. 
Theorem 2.6. There exists exactly one element in 

Proof. In the above proof, if we further assume that u and v are in G^'^°, then we have 
w G A*™°, and hence S[w] > {S[u] + S[v])/2. Thus we must have Dk^n = = Ds = 0. 
From Dn = and = we conclude that u = v. □ 

Besides uniqueness, another particular feature for is that the element in it 
doesn't depend on the value of q. This is due to the fact that Szee equals the constant q 
on A*™° and plays no role in the minimization of In the following we will use z^^ to 
denote this two-component ground state corresponding to magnetization M. 

Remark 2.2. The above proof of uniqueness by the convexity property of Lemma [231 is a 
standard one. Unfortunately, we cannot obtain uniqueness of u g G by the same method, 
at least not in an obvious way. The problem comes from the term 2/3sMq(mi — m i)^ in 
Hs{u). In fact, it is shown in [201 that uniqueness fails for M = q = 0. On the other 
hand, for M = 1, or M = and q > 0, the ground state reduces to a single component 
(cf. Remark IZT]) . and uniqueness can also be obtained by the convexity property. (And 



obviously such one-component ground states are also independent of q.) Except for these 
degenerate situations, however, we do not know how to prove or disprove uniqueness. 
See Section 6.1 for more discussion on the difficulty of proving uniqueness. 
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3 Redistribution and Some Further Properties 



In this section we establish some basic results that will be useful. We'll frequently use 
the method of "mass redistribution", or simply redistribution, introduced in [l20l . For 
convenience we recap the idea below. 

Let /i, /2, fn, gi, 92, dm be nonnegative functions in H^{Vt), where is a domain 
in R^, and ra, m are arbitrary finite numbers. Then we say {gi,..., gm) is a redistribution 
of (/i,...,/n) if gf = for each i = l,...,m, where the coefficients Oy are 

nonnegative constants satisfying Yl^i ^ij = 1 for each j = 1, ...,n. If this is the case, 
it's easily seen that 



i=l j=l 



Moreover, we can prove 



J]|V,.|2< J]|V/,|2. (3.2) 

In the remaining of this paper we shall use the notations Am, G^.g and Eg{M, q) to 
specify the values of M and q explicitly. Redistribution provides a simple and concrete 
way to variate an element in Am into another element, in the same space or in another 
Aa/'- Indeed, from (13.11) . if v is a redistribution of some u g Am, then |v| = |u|, so the 
first constraint A/'[v] = 1 is satisfied automatically, and one needs only to take care of 
the second constraint. Also, the two relations (13.11 ) and (13.21) of redistribution enable us 



to easily compare £[\\\ with £[v]. Precisely, we have £pot[v] = £pot[vi] and Sn[v] = Sn[u] 
from (13.11) . and Skin[^] < ^^fcm[u] from (13.21) . As will be seen, these features make it easy 
to deduce some facts by using redistribution, which might otherwise be harder to obtain 
or need more elaboration. 



3.1 Continuity and monotonicity of Eg{M, q) 

In this subsection we prove that Eg{M,q) is a continuous function, and is increasing 
in each variable. Since the two variables are of quite different natures, we treat them 
separately: Consider Eg{-, q) for fixed q, and consider Eg{M, ■) for fixed M. 



3.1.1 Eg as a function of M 

In the following we fix a g > and consider Eg{-,q). The proof of continuity will rely 
on the monotonicity, and hence we prove the latter first. For this we need the following 
lemma. 

Lemma 3.1. S is bounded on Uo<M<i^M,g- 
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Proof. The assertion is equivalent to say that we can choose for every M e [0, 1] an 
f^"^ G Aa/, such that £[i^] is uniformly bounded in M. This is easy to do. For example, 
let / be any nonnegative function in f] Lyf] such that J f'^ = 1. Then for each 
M G [0, 1], let = 0, We have f^^ g Am and 

£:[f^] = J ||V/|2 + l^/2 + /3„/4 + /3,MY + g}, 

which is bounded by the finite number □ 
Proposition 3.2. Eg{-, q) is strictly increasing on [0, 1]. 

Proof. Let u g Gm,?- We first consider < M < 1. For small 5 > 0, let u(5) be the 
redistribution of u defined by 



2 



Mi(5)2 = (1 



1 



Then u(5) g A(i_2(5)a/- Since u(5) is a redistribution of u, £km[^{^)] < ^A;m[u]. One can 
also check by direct computation that 



SzeeH - £zeeH^)] = qS / K + U^) > 0, 

and 



£s[n]-SsHS)]=f3sS J{ui-u.i)^[2ul + 4uiu.i + 6{ui-u.i)^] > 0. (3.3) 



Moreover, if 5 > 0, strict inequality holds in (13.31) . To see this, for < M < 1, note that 



> (Corollary 12.41) and the fact (mi — u i)^ can not be identically zero. While for 



M = 1, only Ml > 0, and the positivity of (13.31) is obvious. Thus we obtain 



Egiil - 25)M, q) < £[ui6)] < S[u] = Eg{M, q) 

for each small 5 > 0, which shows Eg{-, q) is strictly increasing on (0, 1]. 

It remains to show that Eg{-,q) is strictly increasing at 0. Let {M^} be a sequence 
in (0, 1), Mn 0+, and let g for each n. By Lemma [3Jl, £[u^] is uniformly 

bounded, and hence Lemma [23] implies there is a subsequence {u"^*^'^)} of {u'^} such that 
^n(fc) ^ ^oo weakly in B for some u°° g Aq, and 



Eg{0,q) < S[u°°] < liminf £[u"W] = liminf E^(M„(fc), g) = inf Eg{M,q). 

k^oo k^oo 0<M<1 

Thus Eg{0,q) < Eg{M,q) for every M > 0. To see why strict inequality must hold, 
assume -^^(0, g) = Eg{M, q) for some M > 0, then since Eg{-,q) is strictly increasing on 
(0, 1], we have Eg{M/2, q) < Eg{0, g), a contradiction. □ 
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Proposition 3.3. Eg{-, q) is continuous on [0, 1]. 

Proof. The ideas of proving left continuity and right continuity are different. We first 
prove right continuity. Let u g for some < M < 1. For small 5 > 0, let \y{5) be 
the redistribution of u defined by 

ui{5f' = u\ + 5ul + 5u^_i 
UQ{5f = {I - 5)ul 
u-i{6f = {l-6)u\. 

Let's use Mg to denote M[vi{5)]. Then = M + 5 /(mq + i). Since < M < 1, mq 
and u 1 cannot both vanish, and hence > M for 5 > 0, and — )■ M+ as 5 — > 0+. 
Now since Eg{-,q)h strictly increasing, we have 

Q<Eg{M5,q)-Eg{M,q) (3.4) 

for (5 > 0. On the other hand, since u g GM,q while u(5) need not be in Gm^.q, we have 

Eg{Ms, q) - i?,(M, q) < £[u{5)] - £[m] 

< {SsH6)] - Ss[n]) + i£zeeH6)] - Szee[n]) 

since Skinl^iS)] < Skinl^]- It's easy to check that by letting 6 -> 0+, (13.51) implies 

limsnp{Eg{Ms,q) - Eg{M,q)) < 0. 

This together with (13.41 ) imply the right continuity of Eg{-, q) on [0, 1). 

For left-continuity on (0, 1], we prove by contradiction. Let M e (0, 1]. Assume there 
is a sequence {M^} in (0, 1) such that which Eg{Mn, q) doesn't converge to 

Eg{M, q). By choosing a suitable subsequence, we can assume without loss of generality 
that Eg{M, q) - Eg{Mn, q) > £ for each n, for some e > 0. Now let {u"^} be such that 

G GA/„,g, and let {u"(^)} and g Km be as asserted in Lemma |231 we have 

Eg{M,q) < f [u°°] < liminf^K«] = lim Eg{M^^^k),q) < Eg{M,q)-e, 
a contradiction. □ 
Proposition 13. 3 1 implies the following important approximation result. 

Corollary 3.4. For any M g [0, 1], we can find a sequence {M„} in [0, 1] and a sequence 
{u"}, G GM„,g, such that Mn — )• M, Mn ^ M for each n, and — )■ u°° in Mfor some 



G UM,q. 

Proof. Let {M^} be a sequence in [0, 1] such that Mn M and M„ 7^ M for each n. 
Let {u^} be such that u'^ g GM,,,q- By definition Ar[u'^] = 1 and M[u^] M. Since 
S[u'^] = Eg{Mn,q), by continuity of Eg{-,q) we also have S[u^] Eg{M,q). Thus 
by Lemma [231 {u"^} has a subsequence {u"^'^)} such that u'^^'^^ — u°° strongly in B 
for some g Ga//,^. The sequences {M„(/,)} and {u"^'^)} satisfy the assertion of the 
corollary. □ 
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3.1.2 Eg as a function of q 

Now we consider the function Eg{M, ■) for fixed M e [0, 1]. For u e B+, let's temporarily 
write £[u, q] instead of £[u] to indicate its dependence on q. The proofs of monotonicity 
and continuity of Eg{M, ■) are much easier than those of Eg{-, q) above, and the proof of 
continuity doesn't rely on the monotonicity. Indeed, if qi > q2 > 0, let u e G'M,qi, we 
have 

Eg{M, qi) - Eg{M, ga) > ^[u, qi] - £[u, ga] = (gi - ^2) / + ul,), 



which implies Eg{M, ■) is an increasing function on [0, 00), and is strictly increasing if 

< M < 1. 

On the other hand, for any gi, g2 > and e GM,qi , G (CM,g2' we have 
{qi - q2) J {iu\f + (ul.f) = £[u\ gi] - £[u\ ga] < EgiM, gi) - EgiM, ga), (3.6) 
and 

Eg{M,qi)-Eg{M,q2)<£[u\qi]-£[u',q2] = {qi-q2) J {{ulf + {ul,f) . (3.7) 

From dMl) and (|17]), and the fact / ({u^)^ + (u'l^)^) < 1 for = 1, 2, we find 

\Eg{M,qi)-Eg{M,q2)\ < \qi - q2\, 

and hence Eg{M,-) is continuous. We summarize these results in the following proposi- 
tion. 

Proposition 3.5. For fixed M g [0, 1], Eg{M, ■) is an increasing and continuous function 
on [0, 00). Moreover, it's strictly increasing ifO < M <1. 



Remark 3.1. -^^(0, ■) is not strictly increasing. Indeed, by Proposition (13.71) below, for 
g > 0, u G Go,g satisfies ui = n 1 = 0. Such one-component ground state is unique and 
independent of g (see Remark [Z2)) . Thus -^^(0, •) is a constant function on (0, 00), and 



hence on [0, cx)) by continuity. 

With the continuity of Eg{M, ■) we can prove the analogue of Corollary 13 .41 for fixed 
M and varied g. The proof is the same as that of Corollary I3.4l by changing the roles of 
M and g, and is omitted. 

Corollary 3.6. Fix M g [0, 1]. For any q >0, we can find a sequence {g^} in [0, 00) and 
a sequence {u"}, u" g GiM,q„> such that qn q, q-n ^ qfar each n, and — )■ u°° in B 
for some u°° G GiM,q- 
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3.2 w_i is no larger than ui 

We show in this subsection that u-i < ui for any u e GrM,q- With the aid of this fact we 
will prove u(x) decays exponentially as |a;| oo. 

Proposition 3.7. Let u e Go,g. Ifq > 0, then ui = m_i = 0. 

Proof. Let v be the element in Aq defined by 

2 2 / 2 I 2 \ /o 
Vi = V_i = {Ui + U_i)/2 

2 2 
Wo = Mo • 

Then S[u] - S[v] = (£fci„[u] - SkinW]) + ^^s[u]. Since v is a redistribution of u, Skin[^] - 
^A:m[v] > 0. Also, £s[u] > 0, and hence S[u] - £[v] > 0. But u g Go,g, so we must have 
S[u] - S[v] = 0. Thus actually Skml^] - ^kinW] = ^s[^] = 0, which implies ui = To 
see why they must vanish, note that 



^[u] = y I iVujf + V^(x)|u|2 + /3„|u|^ + q{ul + ul^) 
> j {|V|u||2 + \/(a;)|u|2 + /3,|u|4} =f[(0,|u|,0)]. 



(3.8) 



Again since u g Go,g and (0, |u|,0) G Aq, we must have £[\\\ = £^[(0, |u|, 0)]. Thus 
equality holds in (13.81) . This can happen only if u = (0, |u|, 0) since ^ iVwjp > |V|u|p 



and g > 0. □ 

Remark 3.2. For M = g = 0, we also have ui = m_i but ground states are not unique. 
The case ui = m i = corresponds to only one possibility. See [|20l . Theorem 4.2. 

Proposition 3.8. For every < M < 1 and g > 0, u g GM,q satisfies m_i < ui. 

Proof. Let v be defined by vi = max(Mi, v-i = min(Mi, and vq = uq. Then 
we have £[v] = S[u]. To check this equality, for the kinetic part Skin one can use the 
formula 

for j = 1, -1, and the fact |V|/|p = |V/p a.e. for any / g H^. The equalities of the 
other parts are obvious. Thus Eg{M[v], q) < Eg{M, q), and we have 

M[v]<M (3.9) 

since Eg{-, q) is strictly increasing. On the other hand, it's also obvious by definition that 

vl~v\>u\-u\. (3.10) 



(13.91) and (13.101) imply v\ - t;^ ^ = u\- u"^^, that is v\- u\ = f ^ ^ - ^, of which the 



left-hand side is nonnegative while the right-hand side is nonpositive by definition of v. 
Thus we really have vi = ui and f_i = n i, which means ?i i < ui. □ 
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Corollary 3.9. Let u g Gm^- IfO<M<l,the Lagrange multiplier A in the GP system 
(2.2) is positive. 

Proof, {l.ld) multiplied by n i minus (2.2c) multiplied by ui gives 

2\uiu-i = V ■ (-n_iVui + uiVu-i) + 2/3<,(Mf - u'i-^){ul + 2miu_i). (3.11) 

Thus 

A / uiu-i = /3s {ul-u^i){ul + 2uiu-i). (3.12) 



To justify (13.121) carefully, we can test (13.111) by (fR{\x\), where (fR : [0, oo) [0, 1] is 
defined by 



and then let R go to infinity. 

Now Miu 1 > by Corollary 12.41 and hence J niw i > 0. On the other hand, by 
Proposition 13.81 we have ul - u^_^ > 0, which cannot be identically zero since M > 0. 
Thus we also have / {uf - u'^^){ul + 2mim„i) > 0, and (13.121) implies A > 0. □ 




It's interesting that with this corollary of Proposition 13.81 we can further sharpen the 



assertion of Proposition 1X8 



Proposition 3.10. For < M < 1 and q >0, u e GM,q satisfies m_i < ui. 

Proof. If M = 1, we have m > = For < M < 1, let w = ui- u-i. Then (2.2a) 
minus (2.2c) gives 

Aw + Qw = — A(mi + — fiw, (3.13) 

where Q = -V - 2/3„|u|2 - 2/3^ [2ul + {ui + u^i)^] - g < 0. Since A > and w > 0, 
by subtracting \fi\w from both sides of (13.131) . we obtain Aw + Qw < 0, where Q = 
Q — l/^l < 0. By the strong minimum principle either w > everywhere or w = 0. But 
w = means ui = contradicting to the assumption M > 0. Thus w > 0, which is 
what we want. □ 

Remark 3.3. Recall that we denote the unique element in Gj^^ by z^^, which is inde- 
pendent of q. Assume < M < 1, then we also have z^^ < zf. This is because z*^ is 



just the element in Gm,o from 11201 



With the aid of Proposition 13.81 we now prove the exponential decaying of ground 
states. The approach by using Yukawa potential is exactly taken from Lemma A.5 of 

m. 
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Proposition 3.11. Let u g 'GM,qfor arbitrary < M < 1 and q>0. For any t > 0, there 

exists Mj{t) (j = 1, 0, -1) such that Uj{x) < Mj{t)e^^\''\. 

Proof. (2.2b) can be arranged as 

(-A + t^)uQ = QoUQ, 

where 

Qo = t'^ + ^-V- 2/3„|u|2 - 2/3,(ui - (3.14) 

Thus 

uo{x) = / Yt{x - y)Qo{y)uo{y)dy, 



where Yt{x) = e"*l^'/(47r|2;|) is the fundamental solution of the operator —A + (also 
referred to as the Yukawa potential. See [17], 6.23). By the assumption (Al), Qq < 
outside a bounded set, say B{Ro), the open ball centered at the origin with radius Rq. 
Thus we obtain 



uo{x) < / Ytix - y)Qo{y)uo{y)dy = e ^l^^l / — -Qo{y)uo{y)dy. 

J\y\<Ro Ay\<Ro 

Thus uo{x) < Mo(t)e~*l^^ where 

Mo(t) = sup / ^— ^QQ{y)uo{y)dy. (3.15) 

For uj, j = 1, -1, we similarly have 

(—A + t'^)uj = QjUj — 2[3sU^{uj — u-j) 
from (2.2a) and (2.2c), where 

Qj = t^ + ^i + j\-V - 2/3^|u|2 - 2l3s{u) - ul^) - q. 
Now since m„i < ui, Qi is also negative outside B{Ri) for some radius Ri, and 



ui{x) = J Yt{x - y) [Qi{y)ui{y) - 2f3sUo{yf{ui{y) - u^i{y))]dy 
< / Yt{x - y)Qi{y)ui{y)dy 



< / Yt{x -y)Qi{y)ui{y)dy. 
J\y\<Ri 

As above we conclude that ui{x) < Mi(t)e"*l^l, where Mi(t) is given by (13.151) with all 
the indices replaced by 1. In contrast, the fact m i < ui makes it difficult to apply 
the same argument to m i. Nevertheless, also since m i < mi, at least we can choose 

M_i(t) = Mi(t). □ 
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Later we will consider ground states corresponding to different values of M and g, 
and hence different Lagrange multipliers. The following observation will be useful. 

Lemma 3.12. For n and A in bounded sets, Mj{t) can be chosen to be independent offi 
and A. 

Proof. Take Mo(t) for example. Assume ;U < C for some constant C > 0. From (13.141) , 
Qo<t'^ + fi-V<t'^ + C- V, and hence Rq can be chosen to be independent of /i (and 
A). Then by Holder inequality and the fact J < J'|up = l, (13.151) gives 

^2t(\x\-\x-y\) ^ \ 1/2 



Mo{t) < sup ( / i^Qoivrdy 



<Ro 16vr2|x-2/|2 

/ r 2t{\x\-\x~y\) o 

<^^P( TTT^ ^(^ +C- V{y)fdy 

xeW \J\y\<Ro J-OTT \x-y\ J 

Mi(t) can be estimated similarly, and the assertion for M follows. □ 



4 Redistributional Perturbation in a Fixed Admissible Class 



Let u e GM,q- We have seen it's sometimes useful to construct a "redistributional pertur- 
bation" u{6) of u. In previous examples (namely Proposition 13.21 and Proposition 13.31) . 
the u((5) are so constructed to be in different Am', in order to compare ground states with 
different magnetizations. In this section we investigate the idea, possibly more natural, 
of perturbing u in the same class. Then, since S[u] < S[u{6)] for each 5 > 0, we have 



5=0+ 



>0, (4.1) 



5=0+ 



as long as the derivatives exist. Here denotes right differentiation at 5 = 0. It 

turns out that the existence of such derivatives need some verification. In the following 
we give two examples. Proposition 14.11 and Proposition 14.21 which will be useful in the 
next section. We first introduce some notations. 

• For 5 > 0, we will use D{u{6)) to denote the difference quotient {H{u{6)) - 
H{n))/5. Thus ^^[u(5)]L = 111115^0+ which equals / ^if(u(5)) I , 



5=0+ 



if we can differentiate it under the integral sign. Similarly, Dkin{u{6)), Ds{u{6)) 
and D^ee(u((5)) are understood to be the difference quotients of the indicated parts 
ofH{u{5)). 



For u e GM,q, we write S{ui, uj) = 
we will use the following formula: 



UiVuj-UjVuiY. When computing £)fcj„(u(5)). 



E 



E 



1/2 



3 3 J 



> 0, where a, (j = 1,0,-1) are nonnegative constants. 
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Proposition 4.1. For < M < 1 and q>0, ifz^ is a ground state, i.e. tM e Qu.q, then 



4/5, 



I zfz^,{zf - z^,){TZ^, - zf) > g(l + M) + 



M\2' 



(4.2) 



where r = (/ zl)/{f z\) = (1 + M)/(l - ikf). 



Proof. In this proof we omit the superscript M of for simplicity. Consider the redis- 
tribution u(5) of z defined by 

ux{Sf = {\-b)z\ 
u^iSf = 5zl + t5z\ 
u.i{5f = {l-T5)zl^. 

It's easy to check n{5) e Am for each small 5 > 0. We compute as 
follows: First, 



^fcm(u(5)) = - j |V«o(5)|2 - (5|V^i|2 + t5|V^_i|2) 



t5'(^i,^-i) 



which is independent of 5, and hence 



d 



— £A;m[u(5)] 



T^(2:i,Z_i) 
Z^+Tzi-^ ' 



5=0+ ^1 

Second, for 5 > in a fixed small neighborhood of 0, it's easy to check 



< C\z\ 



for some C independent of 5. Thus it's valid to differentiate under the integral 

sign, which gives 



6=0+ 



^2Ps < {z\ + t£i)(^i - + {z\ - z\){-z\ + t£i) 



/ 



zxz-x{zx - z-x){tz-x - Zx). 



Finally, £^^ee[u(5)] ^ q\{\- 8) \ z\^{\- tS) J 2;^^] , and we have 



^^Zee[u(5)] 



5=0+ 



q{- z\-T \ £1 =-g(l + M). 



The assertion now follows ^£^[u(5)] \^^f^+ > 0. 



□ 
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Proposition 4.2. For < M < 1 and q>0, any u e G^^q satisfies 



q J ul>l3sj ul{ui - u^if + 
Proof. Let u(5) be defined by 



UQ{6 f = (1 - 28)ul 
It's easy to see \y{6) e Am for each small 6 > 0. Now 




(4.3) 



5=0+ 



lim / Dzee{n{6)) = 2q / 

5— >'0+ 



On the other hand, it's not obvious whether we can differentiate £^fcm[u(5)] and £^s[u(5)] 
under the integral signs. We could avoid this problem as follows. Since J D{u{6)) > 0, 



DzeeiuiS)) > - / Dsiui6)) - / Z^fcin(u(5)). 



(4.4) 



Now Dkin{^{5)) < since \i{5) is a redistribution of u. Also, it's easy to check that 
^ifs(u(5)) < for small 5 > 0, and hence Ds{vi{5)) < for small 5 > 0. Thus, after 
taking liminf as 5 — > 0+, we can apply Fatou's lemma to the right-hand side of (14.41) . and 
we obtain 



2q 4> 



5=0+ 



d 

— Hkin{u{6)) 



5=0+ 



It's easy to check that, divided by 2, (14.51) gives (14.31) 



(4.5) 

□ 



Remark 4.1. Now that the terms of the right-hand side of (14.31) are finite, Uq{ui — 



U- 



^i)"^ /{uiu-i) and S{uj,uo)/u'^- (j = 1,-1) are integrable. We can use them to find 



suitable bounds of 



u{6))\ and \^Hkin{u{6))\ independent of small 6 > 0, and 



conclude that £^s[u(5)] and Skinl^i^)] can really be differentiated under the integral sign. 
One might suspect such operations of taking differentiation should be valid for all sim- 
ilar constructions of u{6). This is probably true. However, there are cases of which the 
validity is not easy to see. See Section 6.3 for further discussion. 



There is another point of view on what we did above, which leads us to find (14.31) 
is really an equality but not merely an inequality. We discuss it in the following. At 
any rate, a redistributional perturbation u(5) is a kind of perturbation, and it's natural to 
suspect that the results above could also be obtained from the GP system (2.2), which 
consists of information from general perturbations. Indeed, using chain rule formally we 
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have = £^'[u](u'(0+)), and one would expect (|4.1|) might be a consequence 

of testing (2.2) by u'(0+). This inference is not totally rigorous. Most importantly, we 
are not sure whether u'(0+) is good enough (for example in B) so that £^'[u](u'(0+)) 
makes sense. It turns out that we can really prove equality holds in (14.31) by using the GP 
system, while in our argument the inequality itself plays a critical role. We demonstrate 
this claim in the rest of this section. For this we first give a lemma. 



Lemma 4.3. Let f : 



and g = (91,92,93) ■ 



be continuous. Assume 



\f{x)\ = 0{\x\ ^) as \x\ — )■ 00, and |g| G L^. Then ifV ■ (/g) (as a distribution) is in L^, 
we have J V ■ (/g) = 0. 



Proof. Write V ■ (/g) = h. Then 



h = lim 



/g- n, 



(4.6) 



dB{r) 



where n is the unit outer normal on dB{r). Let's denote | /^^^^^ /g ■ n| by /(r), and our 
goal is to prove limr.^oo/(?^) = 0. By assumption, for |x| large, \ f{x)\ < C\x\^^ for some 
constant C > 0, and hence for r large enough we have 



/(r) < C\r\ 



\g\<C\r\-'(f \g\'Y\dBir)\'/^ = c( [ \gA'^' 

JdB{r) \JdB{r) J \J dB{r) J 



by Holder's inequality, where \dB{r)\ denotes the area of dB{r) and C = C\dB{l)\^/'^. 
As a consequence, we have 



|g|2 ) dr > C-^ / I{rfdr, 

dB{r) J Jo 



(4.7) 



which is finite since |g| e L^. Note that by (14.61) , liiiir^oo ^(^) does exist, and it must be 
zero by (gj]). □ 



It's also convenient to give the following computational result first. The proof is 
straightforward and omitted. 



Lemma 4.4. Assume f,9 e and f,9 > 0, then 

which is regarded as an identity of distributions. 

We can now prove the promised result. 
Theorem 4.5. The inequality (14.31) is an equality. 
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Proof. Let u((5) be as in the proof of Proposition 



O Then 



u'(0+) = {ul/{2ui),-UQ,ul/{2u^i)) . 

From the discussion before Lemma |43] this indicates we do the following: Multiplying 
(2.2a) by Mq/(2ui), multiplying (2.2b) by -uq, and multiplying (2.2c) by Uq/(2m_i). 
Summing these three equations, and after some rearrangement, we obtain 

,ul = P,ul(u, - u^,f ( 2 H- +lTul(^-^). (4.8) 



UlU-l I I — ^ \ Uj Uq 

Now by Lemma 1131 for j = 1, -1 we have 
Note that {ujV {uo/uj)\^ is just S{uj, uq)/u^-, and hence it remains to show 





(4.9) 



m 



V • UQUjV — 

For this we check the conditions in Lemma |431 with f = uq and g = UjV{uo/uj). First, 
from Proposition ITTTl uo{x) = 0{\x\^^). Then, from (14.31) . S{uj,uo)/u'^- is integrable. 



which means \ujV {uo/uj)\ e iP-. Finally, from (l48l) . (l49l) and (1431) . V ■ {uQUjV{uQ/uj)) 
is really in L^, and the proof is completed. □ 

Remark 4.2. To eliminate the unwanted term V ■ {uoUjV{uo/uj)), in the proof above 
we use the inequality (14.31) at two places: to guarantee that \ujV {uo/uj)\ e L^, and to 
guarantee V ■ {uQUjV{uQ/uj)) e L^. It looks somewhat pedantic, but seems unavoidable. 
For example, from (14.81) alone, we do not even know if Uq{ui - m_i)^(2 + Uq/{uiu^i)) 
is in L^. We remark that similar problems happen to other constructions of u(5) which 
lead to equalities. Thus in our context the inequalities obtained from redistribution are 
not consequences of the GP system. This declaration however may be overthrown if we 
can prove some comparison results of the asymptotic behaviors of the three components. 
See Section 6.3 for discussion. 

Remark 4.3. We'd like to do the same thing for (14.21) . However, note that since z^^ is 
independent of q, it's impossible that (14.21) be an equality for varied q. Indeed, following 
the idea of proving Theorem 1431 we get a trouble at the very beginning: by letting u(5) 
be defined as in the proof of Proposition 14. 1[ we have u'(0+) = +oo, -tz^\/2), 

which suggests we multiply (2.2b) for z*^ (which is the trivial equation = 0) by infinity. 
This problem can be avoided if there is a sequence of ground states (corresponding 
to different values of M or g or both) such that Un > 0, and Uq tends to zero. The details 
will be given in Section 5.2. 
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5 Bifurcation Between 2C and 3C Ground States 



In this section we study the bifurcation between the 2C regime and the 3C regime of 
ground states for (M, q) e (0, 1) x (0, oo). (We recommend the reader Figure 5 of [|T9l 
for a clear illustration.) In [20], we proved that for < M < 1, u g Gm,o implies uq = 0, 
i.e. u = z^^ . According to numerical results, if q is not too large, ground state remains 
2C, while for q greater than some critical value no emerges. We prove this observation 
in Section 5.1. Some characterizations of the bifurcation points are given in Section 5.2. 



5.1 Existence of bifurcation point 

The observation we are going to prove is summarized as the following theorem. 

Theorem 5.1. For < M < 1, there is a qc{M) > such that for q > qdM) (resp. 
q < qc{M)), u G GM,q implies uq > (resp. u = z^). 



The proof is separated into several parts. Roughly speaking, our idea is using (|4.2I) 
to disprove uq = 0, and using (14.31) to disprove uq > 0. In the following we fix an 

M G (0,1). 

Claim 1 For q large enough, if u g GM,q, we have uo> 0. 

Proof. Assume z^^ g G'M,q, then z^^ satisfies (14.21) . Since z^^ is independent of q, it 
indeed gives an upper bound for g, and hence the assertion. □ 

Claim 2 Assume for some q there exists u g (Gm,^ with uq > 0, then for every q' > q, 
V G GM,q' satisfies fo > 0. 

Proof. Let's use S[-,q] instead of £[■] to specify the value of q. Since u g £[u, q] < 

£[z^'^,q]. Thus, by the assumption uq > 0, we have S[u,q'] < S[z^'^,q'] for any q' > q. 
Hence z^ ^ Ga/,^' for q' > 0, which is what we want to show. □ 

Now define 



qc{M) = inf |g | for any q' > q,v e GM,q' have no > } . 

From Claim 1, qc{M) < oo. By definition for any q > qc{M) and v g Gm,^, we have 
no > 0. Moreover, since qc{M) is the infimum of all numbers satisfying this property. 
Claim 2 implies qc{M) also satisfies the following property: For any < g < qc{M) and 
V G (Gm,(j» '^0 = 0. To complete the proof of Theorem 15. 11 it remains to show qc{M) > 0. 
This is the most difficult part, for which we need to prove the following result. 

Claim 3 There exists g > such that u g (GM,g implies u = z^. 
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To disprove the presence of uq, we shall use (14.31) . In fact we will only use q J Uq > 



2(3s J ul{ui - M-i)^, from which it's easy to see = if g = 0. This is the argument 
used in [20]. For g > 0, however, whether = is not so obvious. Our proof of Claim 
3 is based on the following lemma. 

Lemma 5.2. Given q > 0, if qn ^ q and g GiM,q„ is such that converges to some 
u e GiM,q in B, then the following assertions hold. 

(a) There exists R> independent ofn such that 

[<? < [ (5.1) 

Jb{r) 

(b) u" — )■ u uniformly. 

We first prove Claim 3 by these two lemmas. 

Proof of Claim 3. By Corollary 13.61 there exist a sequence qn -> 0+ such that gn 7^ 
for each n, and a sequence u'* g GM,q„ such that u"^ -> z^^ in B. Let R be as asserted 
in Assertion (a) of Lemma 15^ and let k = inf5(^)(zf^ — z^\). Note that A; > by 
Remark 13.31 Now by Assertion (b) of Lemma 15.21 u'* — z^^ uniformly, and hence 
-u'^i) > k/2 on B{R) for large n. From this fact and Assertion (a) we obtain 

Jb{r) 

> ^ / (5.2) 



<B{R) 



A-2 , 



for n large enough. On the other hand, for any n, (14.31) implies 



qn J {u^r > 2/3. J (n^)^K - • (5.3) 

Since gn 0+, (15.21) and (15.31) implies Uq must vanish for large n, which completes the 
proof. □ 

Now we prove Lemma 15^ 

Proof of Lemma \5?2\ We first prove Assertion (b). The idea is that, if the GP system 
(2.2) for u"^ tends to that for u, then uniform convergence can be obtained by the global 
boundedness result for elliptic operators. For this purpose, we need to show that the 
Lagrange multipliers for u'*, denoted by and A^, converge to those for u, denoted by 
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yu and A as before. This can be done as follows. Multiply (2.2a) by ui and multiply (2.2c) 
by u-i, and integrate, we obtain (/x + jX) J = Fj{u, q) for j = 1,-1, where 

F,in, ^ I { I V,,f . + 2A.|uP„| + - ,._,) + u^iu^ - + ■ 

Solve for /x and A we obtain 



/z = [Fi(u, q)/{J ul) + F„i(u, g)/ / i] /2 
A = [Fi(u, ul) - F_i(u, g)/ / ul,] /2. 



(5.4) 



/in and An can also be expressed by the above formulas, with u replaced by and 
q replaced by qn- Note that since we assume < Af < 1, each J (wp^ and J uj are 
nonzero for j = 1,-1. As a consequence, fin ^ and An -> A follow the fact u"^ u in 
B. 

Now let = - Uj. Subtract (2.2a) for u from (2.2a) for u'^, we obtain 

Av'l-Vv'l = Pn + S{u'')-S{u), (5.5) 



where 



Pn = -(/in + An - qn)Ul + (yU + A - q)ui, 

S{u) = 2/3n|u|\i + 2(3s [ul{ui - + ui{ul - uii)] . 



Apply global boundedness theorem for elliptic operators (see e.g. 11101 , Theorem 8.16) 
to (15.51) . we have 

sup Kl < sup Kl + C\\Pn + Siu"") - S{u)\\l2, (5.6) 

B{r) dB{r) 

where C > depends only on the radius r and sup5(j.) V . Now since /i, An A, 
qn q, and -> u in B, we see Pn in L^, and also S'(u"') - S'(u) in 
since is continuously embedded in L^. On the other hand, also since fin ^ and 
An A, by Lemma [37121 we can find Mj{t) independent of n such that Uj and each 
are bounded above by Mj(t)e^*^l. In particular, for any e > 0, there exists r > such 
that 

supK(x)|<£ (5.7) 

|a::|>r 

for all n. Fix this r in (15.61) and let n — )■ cx), we obtain 

limsup ( sup |'yi*(a;)|J < e. (5.8) 
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From (15.71) and (15.81) we have supRs | < Is for large n. Since £ > is arbitrary, 
we conclude that f ^ — )• uniformly. Similarly we have and v^-^ converge to zero 
uniformly, which completes the proof of Assertion (b). 

Then for Assertion (a). As mentioned we have /i and hence ii^ forms a bounded 
sequence, say iin<C for some constant C > 0. Now multiply (2.2b) for u"^ by Mq, and 
then integrate, we obtain 

\^n^{ulf = I {|V«S|2 + VK)2 + 2/3,|u-|2(«-)2 + 2/3,K)2K-«-,)2}, 



which implies 



I V{u^f <f^nj <C j (5.9) 



On the other hand, by the assumption (Al), there exists > such that V{x) > 2C for 
|x| > R, and hence 

Vk)2>/" V{u^f>2c[ (5.10) 
Jb{ry jb{ry 



From (15.91 ) and (15.101) , we obtain j' [uqY > 2 JB(R)ci'^o)'^^ which is equivalent to (15.11) . 



□ 

Now we have completed the proof of Theorem 15.11 We remark that, however, it 
doesn't provide a good description of qc{M), even in a qualitative sense. Notably, we 
don't know why qc should be a continuous increasing function of M, as is quite apparent 
from the numerical results. Nevertheless, one fact that is not quite clear numerically can 
be settled by our method. That is, as M — )■ 1~, whether qc{M) is tending to infinity or 
some finite number. The following complement of Theorem 15. 11 says it's the latter that 
is the case. 

Theorem 5.3. There is a q > such that if q > q, u e GM,q satisfies uq > for every 
M e (0, 1). 



Proof Let's denote the left-hand side of (l42h by L(M), i.e 



L(M) = 4/3, j z^z^.izf' - z^,)irz'{ - zf'). 
It suffices to show that L(M) is uniformly bounded for M e (0, 1). Note that 

L(M) < 4/3, J zfz^, ■ (^f ) ■ (rz^,) = A^^r j {z^f{z^,f 
< 4/3.r||.nii. / {^""i? = 2/3.(1 + M)\\z^\\l^. 
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Thus it suffices to show that 11 uniformly bounded. However, since z*' is the 

unique element in Gj\/o, from Corollary 13.41 the map M z^^ from [0, 1] into B is 
continuous, and we can prove that M is also continuous from [0, 1] into by 

imitating the proof of Assertion (b) of Lemma 15^ which completes the proof. □ 

Remark 5.1. It might be a little surprising that, by the same argument, we have trouble 
to conclude that M h- z^\ is also continuous from [0, 1] into L°°. Indeed, the problem 
only occurs at M = 1, where z^{ is equal to zero. See Section 6.2 for discussion. 



5.2 Some characterizations of qc{M) 

In the following we also fix an M e (0, 1). Although our statement of Theorem 15. II is 
a qualitative one, our proof does provide some quantitative information. For example, 
as the proof of Claim 1 says, (14.21) gives an upper bound of the value q for which z^ e 
GM,q, which means we have an upper bound of qc{M) in terms of z*^. Precisely, we 
have QciM) < U{M), where 



U{M) = ^— { 4/3, / - z'1\){tz'1\ - zf') - 



1 + M\"J ' ' ' J {z^Y + t{z^^Y\' 

Similarly, (14.31) implies qc{M) > inf -Ffu], where 

^3sjul{ul-u^lf{2 + ul/{ulU^l)) + jY.^=l-lS{uJ,UQ)/u] 

J K 

and the infimum is taken over all u g Ug>oGM,g such that uq > 0. Indeed, from Theorem 
14.51 this inf J^[u] is equal to qc{M). To see this, for q > q^M) let be an element in 
G^f g, then > 0. By Theorem |431 we have q = J'[u'^], and hence 

gc(M) = lim q= lim J^[u^] = inf J^[u]. (5.11) 

Of course this characterization is not of much use as to compute qc{M). In contrast, to 
obtain the upper bound U{M), one needs only to find the 2C ground state z^^, which 
really gives a great reduction in computation cost. To conclude this section, we explain 
that we can also modify U{M) to a characterization of qc{M) as long as there is a se- 
quence u"^ G GM,q„, where g„ qc{M)^, such that u'^ z*^ in B. Before doing so, 
we remark that the existence of such sequence is left open in this paper, since we do not 
show that z^^ is the unique element in GM,q^{M)- See the discussion on uniqueness in 
Section 6.1. 

Now for q > qc{M) and u g GM,q, consider u(5) defined by 




(1-5)^2 



1 

(1 - a5)u\ . 
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Here a = a{u) = {J uf) / (J u^^). This u{6) can be regarded as the same as that given in 
the proof of Lemma WT\ except for uq > 0. In particular note that a-^r = (l + M)/(l- 
M) as J" Mq 0. One can check that such defined u{S) is in A^^. Follow the idea of 
proving Theorem 14.5 1 we will get the following equality. 

Theorem 5.4. Let M e (0, 1) and q > qc{M). u g Gm,^ satisfies 

4/3s j - M_i)((TM_i - Ml) + 2/3s j ul{ui-u-i){au-i-ui) 

S{uQ,ui) + aS{uQ,u^i) 



Now if we have qc{M)^ and u"^ g GM,q„ such that z^^ in B, by the above 
theorem we get the following characterization of qc{M): 



qc{M) 



1 + M 



S{u^Q,ul) + anS{u^Q,u\] 



— lim 



(5.12) 



n\2 



where cr„ = (J" (m^)^)/(J" (m"^)^). It remains to show that (15.121) gives the upper bound 
U{M). To see why (15.121) gives the upper bound U{M), let f = Vmq/mo (for general 
u G B with Uq > 0), and we have 



S{uo,ui) + ^^(mo,^-!) 



IVmi— Mifl + cr |VM-1 — M_if I 



u, 







+ (TM_i)|fp - 2(miVmi + (tm_iVm_i) • f + IVmiP + Cr|VM_l| 

2 



> 



/ 2 I 2 N 



aS{ui,U-i) 



f - 



MlVui + CTM-iVM-I 



au 



aS{ui, M_i) 



(5.13) 



Thus, by choosing a subsequence u"(^) of u'* so that u'^^'^) — )• z^^ and Vu"^'^) — )■ Vz^ 
almost everywhere, we can apply Fatou's lemma to obtain 



lim 



n{k) n{k) 







U 



rtiKiN . r,/ n(k) n(k).. 

1 ) +(^nik)S{uQ' ',U_\') 



[U. 



n{k)^2 



> 



TSiz^,Z^^ 



yM\2 



+ nz 



M\2' 



By substituting this inequality into (15.121) we find qdM) < U{M). 



6 Discussions and Open Problems 

In this section we discuss some natural questions arising from this paper. They are 
categorized into three subsections. 
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6.1 Uniqueness 

Uniqueness is a standard and prominent problem to be settled in variational problems. 
In this paper, although it's not essential, from time to time our presentation was plagued 
by the lack of it. For example, assume we have uniqueness, then we can simply use a 
symbol u^^'^ to denote the element in GM,q, and the wordy statements of Corollary |3]4 
and Corollary 13 .61 simply say that (M, q) u^^'^ is a continuous map from [0, 1] x [0, oo) 
into B. Nevertheless, we have mentioned in Remark 1X2] that our energy functional S 
doesn't have the suitable convexity property due to the term 2(3sul{ui - appearing 
in Hs{u). As to remedy this difficulty, there are two natural ideas: 

(a) Although £ is not convex on B, it might be convex on a fixed A^, which is suffi- 
cient to prove uniqueness. 

(b) In this paper there is no assumption on the magnitude of f3s, while in practical 
spin-1 BECs (see references given in the introduction) it's very small compared to 
/3n, and hence Sg contribute to a rather insignificant amount of the whole energy. If 
we are willing to take this into consideration, maybe the convexity of other parts 
will outweigh the nonconvexity of £s. 

We are here to show that both ideas do not work. A counterexample is as follows. 
Let /, g, h be any three nonnegative functions in n Ly n L'^ such that 1) /, g and h 
are supported on disjoint sets, 2) J (/^ + g'^ + h'^) = 1, and 3) J g^ = J h'^ > and 
/ (/^ - 5'^) = M. Then let u = (/, g, h) and v = (/, h, g). We have u, v, w g , where 
Wj = {u'j + Vj)^/'^ as in the proof of Lemma 1231 Let = supp(/), Qg = supp(5f) and 
Qh = supp(/i), then it's easy to check that 




Thus, no matter how small (3s is, £ doesn't have the desired convexity property on Am- 
Of course, the u and v above are far from being ground states, especially due to the 
assumption that the supports of the components are disjoint. We can go on to suspect 
8 might satisfy ^(^[u] + £[v]) - £[w] > as long as u and v are "similar to" ground 
states. Anyway, uniqueness for our model, if holds, can not be obtained from the usual 
convexity argument. On the other hand, it's also not quite clear whether uniqueness 
holds from numerical simulations. The trickiest part lies on the bifurcation point. To 
have a better understanding of the problem, remember that the "nonuniqueness" point 
(M, q) = (0, 0) connects two boundary regimes which sharply contrast to each other: For 
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M > and g = 0, u g Gm,o has uq = 0, while for g > and M = 0, uq is the only 
nonvanishing component of u g Go,q. It's observed in numerical simulations that such 
sharp contrast also occurs along the curve of bifurcation points {M,qc{M)), so sharp that 
one is not easy to tell whether uq shrinks to zero rapidly as q decreases to qc{M), or 
indeed there are both 2C ground state and 3C ground state at qc{M). In ffT9ll , the latter is 
claimed to be the case. However, in other simulations by using numerical continuation 
method (from not published private discussion. See [7] for related study of bifurcation 
with respect to and 15s), one can really track the changes of ground state from 3C 
profile to 2C profile as q decreases from a large number to zero, and it seems ground 
state is always unique (for (M, q) ^ (0, 0)). 

6.2 Uniform variations of ground states at boundary regimes 

We have stated the bifurcation phenomenon in terms of varying q and fixed M. This 
choice is physically natural as the value of q can be tuned by modifying the applied 
magnetic field. From a mathematical point of view, however, we'd like to remark that 
this choice is in fact made intentionally. Somewhat unexpectedly at first sight, there are 
two difficulties to imitate the proof of Theorem 15. II if we consider the bifurcation with 
varying M and fixed q. The first one is that we lack an analogue of Claim 2 in Section 
5.1. Precisely, we do not know how to prove that if there is u g Gm,^ satisfying uq > 0, 
then every v g Gm',? with M' < M must have vq > 0, or equivalently z^^ g GrM,q 
implies g Gm',^ for M' > M. As a consequence, we can't conclude that there 
exists a number Mc{q) which definitely separates the 2C regime and the 3C regime. 
The second problem, more fundamental, is that we are not sure whether the Lagrange 
multipliers will converge as M tends to 1~ or 0+. Note that in either case J ^ — )• for 
u G GA/,g, and we can not use the formula (15.4b directly. As a consequence, we can't 
obtain uniform convergence when M — )■ 1"^ or 0+ as in Lemma \5?2\ Despite of this, 
we remark that in either situation it's known that the component which is not tending to 
zero does converge uniformly. For example, let — > 1 and g 'GsM„,q converges in 
B to the unique element in Gi,g, which we denote by u°° = {uf", 0, 0), then we also have 
u\ uniformly. This is because anyway /i^ + An converges, and (2.2a) for u"^ tends 

to (2.2a) for u°°. What left open is whether Uq and converge to zero uniformly. This 
lack of uniform convergence (of precisely) then prevents us from imitating the proof 
of Claim 3 to conclude that Uq = for large n. Similarly, when M — )■ 0+, we only know 
UQ converges uniformly but not for m and (Of course, this is sufficient to conclude 
that Mo > when M is close to zero.) As we have mentioned in the remark after Theorem 
15.31 such problem also occurs for z^^ when M 1 , where z^-^ converges to zero in B, 
and we don't know if it converges uniformly. 
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6.3 Asymptotic behaviors at infinity 

We are not sure whether ^£^[u(5)] exists for some constructions of redistributional 
perturbation u{6). An example is given by 

Mi(5)2 = {1 - 6)ul + a6ul^ 

M_i(5)2 = 6ul + {1 - a6)uli , 

where assume u e GM,q is such that Uj > for each j, and cr = (J" u\)/{j u^_i). It's 
easy to see u{6) e Am- Thus we have D{u{6)) = Dkmi'^i^)) + Ds{u{6)) > 0, but we are 
not able to deduce from it that ■^Skin[^{^)]\g_Q+ + is^s[u{6)]\^_^^ > since we do not 
know if the derivatives exist. Indeed, if we differentiate them formally under the integral 
signs, we obtain 

aS(u^ + < 2ft. / (ul - ul,)(aul, - „f ) f ^ + 2) . (6.1) 

while neither the left-hand side nor the right-hand side are obviously finite. Note that for 



the left-hand side of (16.11) . we know the term S{ui,u-i)/u\ = |Vm_i - Vwip 
is integrable since m i < mi, and it's ^(mi, m_i)/m^^ that causes trouble. The problem 
here is very similar to that mentioned in Remark 521 Roughly speaking, they are all due 
to the fact that we do not have a comparison of the asymptotic behaviors of different 
components of ground states. We remark that some numerical experiments show that 
uq{x) < u^i{x) < ui{x) as |x| large. In fact it looks like uq/u^i and n_i/Mi 
as I a; I oo. If this can be proved, then the right-hand side of (16.11 ) is finite, and we 
can justify the differentiation above. Also, one can see that the problem mentioned in 
Remark 14.21 disappears, and Theorem 14.51 can be obtained from the GP system (2.2) 
without using Proposition 14.21 



Appendix 

We prove Lemma 1231 in this appendix. For convenience we restate the assertion below. 



Lemma. Let {u^} be a sequence in 1B+. Suppose A^[u"^] 1, A^[u"] M, and £[n^] 
is uniformly bounded in n, then {u"^} has a subsequence {u"^^^)}^^ converging weakly 
to some u°° G A, which satisfies S{\i°°] < \\m'm.ik^oo £[vi^^^\ If we assume further that 
£[vJ^] Eg, then u°° G G, and u''(^) u°° in the norm o/B. 

Proof Since the five parts Skin, Spot, etc. are nonnegative functional on B, the uniform 
boundedness of S[u^] implies that of the five parts, and hence {u"} is a bounded se- 
quence in B. Thus {u"} has a weakly convergent subsequence {u"(^)}^^ in B, of which 
we denote the weak limit by u°°. 
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We first prove e A. Since {u'^'} is a bounded sequence, J \/|u"(^) p < C for some 
C > independent of n. By the assumption (Al), for any e > 0, there exists Re > such 
that V{x) > C/e for |a;| > Re. Thus we have 

C> IV\U<^\^> [ V^|u"W|2>£ f |^n(fc)|2^ 

J Jb{r^y ^ Jb{r^y 

and hence j^^^ y |u"^(^)p < e for each k. In particular j^^^ )c(^j^^^)^ < ^ for each k and 
j = 1, 0, —1. The weak convergence of u"'^'^) in (//^(M"^))^ impHes u^^^^ weakly in 

l2(M"^), and u^^^'* strongly in L^{B{Re)). Hence we have 



j (wf )2 <limmf j ( 



< limsup / {u^^'^^)'^ 



limsupf / {uf^f+ / (uf^)' 



(6.2) 



<e+ {ufY <e+ {uf)\ 

Jb{r,) j 

where the first inequality is due to the weak lower semi-continuity of a norm. Since 
e > is arbitrary, the first and the second inequalities of (16.21 ) must be equalities, which 
imply 

^lim I {uf^f = I {uf f (6.3) 

for j = 1, 0, —1, and hence u°° g A. 

Now since u°° is the weak limit of u'^(^), S[u°°] < liminf £^[u"^('^)] is a consequence 



of standard weak lower semi-continuity theorem. See e.g. Theorem 1.6 of 11271 . Indeed, 
by the same theorem we have 



iVwf P < liminf 1 


' IV7 ni>^)l 


■' k^oo / 





Viuff < liminf / V{uf^^ f, 

fc— !>0O J 

and 

for every continuous function / : [0, oo). In particular every parts of £ satisfies 

such lower- semicontinuity inequality. We claim that these liminf s are all limits and the 
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inequalities are all equalities provided £^[u'^('^)] Eg. This is easily seen by assuming 
otherwise. For example assume f |Vm°°P < limsup f iVwJ^'^^p for some j. Then we 
obtain 

Eg= lim ^[u"^'^)] >S[u°°], 

fc— >-00 



contradicting to the fact u°° e A. Thus the claim is true. This fact together with (16.31 ) 
imply — )■ Now since is reflexive, u"'^'^) ^ u°° weakly in 

and — implies u"^'^) — )■ u°° strongly in H^. Similarly we can prove 

^n(fc) _^ ^oo 2.2, and in L^. □ 
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